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ABSTRACT 

We use N-body simulations to survey the response of embedded star clusters 
to the dispersal of their parent molecular cloud. The final stages of the clusters 
can be divided into three classes: the cluster (i) is destroyed, (ii) has a loose 
structure, and (iii) has a compact core. We are interested in three of the governing 
parameters of the parent cloud: (i) the mass, (ii) the size, and (iii) the dispersing 
rate. It is known that the final stage of the cluster is well correlated with the star 
formation efficiency (SFE) for systems with the same cluster and cloud profile. 
We deem that the SFE alone is not enough to address systems with clouds of 
different sizes. Our result shows that the initial cluster-cloud mass ratio at a 
certain Lagrangian radius, and the initial kinetic energy are better indicators for 
the survivability of embedded clusters. 

Subject headings: methods: A^-body simulations - open clusters and associations: 
general - galaxies: star clusters 



1. Introduction 

Stars are the fundamental units to build up a galaxy and are mostly not born alone 
but in groups. Stellar groups, clusters, or associations are formed inside molecular clouds. 
They are embedded and not optically seen until they get rid of the residual material after 
star formation. An embedded cluster is not always understood clearly because of observa- 
tional constraints. The near-infrared all sky s urvey, Two Micron A ll Sky Survey, has given 



astronomers a chance to look into the clouds. iLada fc Ladal (120031 ) claimed that the num- 
ber of embedded clusters declines rapidly with age and the infant mortality of clusters (i.e., 
disruption not long after birth) is more than 90% in our Galaxy. 
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Galactic tidal forces, close encounters with giant molecular clouds fiGieles et al.ll2006l ). 
shock heating, and mass loss of the syste m by massive member s tars, such as UV radiation, 
stellar winds, and supernova explosions ( Boily fc Kroupa 2003a| jbl). are possible disruption 
mechanisms for stellar groups, whether they are embedded or not. Nonetheless, most of these 
mechanisms (such as the first three aforementioned mechanisms) have a destruction timescale 
longer than the u pper limit of the li fetime of rnolecu l ar clouds, which is about a few to a 



few t ens of Myr (IBlitz &: Shul Il980l : Elmegreenl I2OOOI : iHartmann et all I2OOII : iBonnell et al. 



20061 ). CO observati ons in our Galaxy sug gested that the lifetime of molecular clouds is of 



the order of 10 Myr (ILeisawitz et al.lll989l ). Timescale estimation of photoevaporation and 



statistics on the expected numbers of stars per cloud sho wed that giant molecular clouds 



of mass IO^Mq are expected to survive for about 30 Myr (j Williams fc McKed 119971 ). The 



most promising disruption mechanism for embedded clusters seems to be the dispersion of 



evolution (see, e.g.. 


Tutukov 


1978; 


Lada et al. 


1984; 


Goodwin 


1997; 


Boily & Kroupa 


2003b 


Baumeardt & Kroupa 


2007; 


Bastian & Goodwin 


2006; 


Goodwin & Bastian 


2006 


)• 



In a large set of simulations, iBaumgardt fc Kroupal (120071 ) studied the dispersal of the 
residual gas by decreasing the mass with different star formation efficiency (SFE), and in 
different tidal fields. They concluded that the clusters had to form with SFE > 30% in 
order to survive gas expulsion, and the external tidal fields have significant infiuences only i f 
the ratio of half mass radius to tidal radius is larger than 0.05. iGoodwin fc BastianI (120061 ) 
and iBastian fc GoodwinI (120061 ) addressed a similar problem and found that the embedded 
clusters would be destroyed within a few tens of Myr if the "effective star formation efficiency" 
(hereafter eSFE) < 30%. 



Chen fc Ko I (120081 ) studied the behavior of embedded clusters when the parent cloud 
is dispersing (i.e., the size increases but the total mass remains constant, thus the cloud 
becomes more and more diluted). From a large survey of simulations we found that the 
final structure can be classified, according to the expansion ratio of r^^ (the 45% Lagrangian 
radius of the cluster), into three groups: (i) destroyed, (ii) loose structure, and (iii) compact 
core. Empirically, the expansion ratio is related to the cluster-cloud mass ratio (which will 
be defined later in this paper). 

In this work, we address the problem of what determines the fate of the embedded 
cluster. We consider that the cluster and the parent cloud could have different initial density 
profiles. Hence we examine three parameters of the parent cloud independently: (i) mass, 
(n) size, and (in) dispersing rate of the cloud. Clusters with and without mass functions 
are considered. The paper is organized as follows. In Section 2, we describe the model and 
simulations. In Secttion 3, we discuss the applicability of SFE as an indicator for describing 
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the behavior of the cluster, and then introduce other better indicators for more general initial 
conditions. A summary and some remarks are provided in Section 4. 



Model 



To study the survivability of embedded star clusters, we consider an idealized model. 
We put a spherical star cluster at the center of an expanding spherical molecular cloud, and 
then examine the su bsequent behavio r numerically. In our simulations, the distribution of 
stars in the cluster is iPlummerl (jlQllI ) with a length scale of 0. 6 pc (note that the half mass 
radius is ~ 0.8 pc). The cloud is represented by an external iPlummerl (jlQllI ) potential in 
which the length scale increases with tii ne to represen t the dispersion of a cloud. We adopt 
the A^-body simulation code NB0DY2 (jAarsethll200ll ) to study the evolution of the clusters. 



Similar systematic sirnulation s have been performed by , e.g., iGever fc BurkertI (120011 ) 
and Baumgardt &: Kroupa ( 2007 ). Geyer fc BurkertI (2001) used the King ( 1966) model 



to des cribed both the cloud and gas, and iBaumgardt fc Kroupal (l2007l ) used the iPlummer 
(llQllI ) model. Both papers considered equal-mass stars and represented the cloud dispersion 
as a continuous reduction of cloud mass. In contrast, we mimic cloud dispersion as the 
expansion of the cloud and the cloud mass does not change. 



2.1. Initial condition 



We assume that each embedded cluster contains 2500 stars with a total mass of 2500 
Mq. The physical distribution is Plummer and the length scale, ac, is about 0.6 pc. The 
velocity dist ribution depend s on the potential of the parent cloud. The initial mass function 
(IMF) with ISalpeterl (119551 ) slope from 0.3 to 30 Mq is also considered, and there is no 
primordial mass segregation. All clusters, with the parent clouds, are required to be in virial 
equilibrium prior to cloud dispersion. 

The SFE (and more precisely, the total SFE) is defined as 



Mr 



V 



Mr + Mb 



where is the mass of the cluster and Mh is the mass of the parent cloud. In our runs M^ 
ranges from 0.5 to 19 M^, which gives t] from 67% to 5%. 



GoodwinI (120081 ) also discussed cases of a cluster and cloud with different spatial distri- 



butions but the cluster is in virial equilibrium initially. This is very similar to how we set up 
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our initial conditions. Note that the eSFE of iGoodwinI (120081 ) is related to the length scale 
of the cloud. More discussion on eSFE will come later. 



2.2. Cloud dispersion 

The cloud is represented by an external Plummer potential, 



(2) 



where G is the gravitational constant, r is the distance from the cloud center (which is also 
the star cluster center), and a is the length scale of the potential. To mimic expansion, we 
let a increases with time, 

a = aoe*/*^ , (3) 

where Cq is the initial length scale and te is the dispersing e-fold timescale. A molecular cloud 
has a relatively short lifetime ranging from a few to a few tens of Myr. We try four different e- 
fold timescales: from 0.33 to 3.3 Myr. The potential becomes ineffective in no more than two 
e-fold times (largest dispersing timescale) to 10 e-fold times (smallest dispersing timescale), 
which corresponds to 7-3 Myr. As molecular clouds are not expected to last longer than 30 
Myr, we stop our simulations at 30 Myr after the cloud starts dispersing. In any case, the 
influence of the gas removal is long gone well before the end of the simulations. 



3. Results and discussion 

We assume that the embedded clusters are bound and in virial equilibrium initially. 
When the cloud disperses, the cluster starts to expand. It may be destroyed or may shrink 
back later by self gravity. 



3.1. Cluster and cloud have the same initial density profile 

In this subsection, we present some cases for clouds having the same Plummer length 
scale as the clusters do, 0.6 pc. In these cases, SFE is a good indicator for the subsequent 
behavior of the clusters. Figure [T] shows the half mass radius (r/i^) evolution in different 
SFEs with a cloud dispersing (or gas removal) timescale of 1.1 Myr. rhm increases almost 
linearly when SFE is less than 20%. Clusters are destroyed in these cases. For cases with 
SFE from 25% to 50%, rhm expands slower and shrinks back (even becomes stable) later. 
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In these runs, clusters remain intact but with different concentrations. The expansion of a 
cluster by gas remov al and then shrinking (colla pse) to equilibrium later w a s also reported 
in the simulations by iGoodwin fc BastianI (|2006[) and iBaumgardt fc Kroupal (120071 ) , and has 
also been compared with observations by iBastian et al.l (120081 ) . 




t [Myr] 



Fig. 1. — Evolution of half-mass radius as a function of time. Here, the dispersing timescale 
te is 1.1 Myr, SFE is from 5% to 50%. 



3.1.1. Different dispersing timescales 

The behavior is expected to be different if we change the cloud dispersing timescale (or 
gas removal timescale). Figure [2] shows the evolution of SFE = 25% with four dispersing 
timescales. increases almost linearly when is 0.33 Myr, and increases slower when te 
is 0.65 Myr. For t^ equal to 1.1 and 3.3 Myr, r^m increases much more slowly and changes 
little after some time. 



3.1.2. Expansion ratio of o,nd hound mass fraction 



After dispersing, all clusters expand. To compare the structure before and after, we 
define the expansion ratio of the half-mass radius rhm as 



rfern (final) 
r/jm (initial) ' 



(4) 
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Fig. 2. — Evolution of half-mass radius as a function of time, in terms of Myr. The dispersing 
timescale here is from 0.33 to 3.3 Myr. The SFE is 25%. 

Here final means 30 Myr. 

Figure [3] represents the results for the expansion ratio of r/^^ at the end of the sim- 
ulations, 30 Myr, as a function of SFE in different dispersing timescales. Different dis- 
persing timescales produce different expansion ratios for the same SFE. For each dispersing 
timescale, there is a simple relation between rj and e. Each relation has two branches jointed 
by a turnover point. Specifically, the turnover points for = 3.3, 1.1, 0.65, 0.33 Myr are at 
rj ^ 0.18,0.25,0.33,0.50, respectively. The steeper branch corresponds to destroyed clusters 
and the flatter branch to survived clusters. Hence the SFE, 77, is a good indicator for the 
survivability of embedded clusters. 

Figure H] depicts the bound mass fraction as a function of SFE in different dispersing 
timescales (or gas removal timescale). Both Figures [3] and H] nicely demonstrate the expected 
result: a cluster might survive in a slowly dispersing cloud but might be destroyed in a fast 
dispersing one. The depend ence of the bound mas s fraction on SFE is als o reported by 



Baumgardt fc Kroupa] (120071 ). Based on this relation. iParmentier et al.l (120081 ) discussed the 



shape of the initial cluster mass function. 
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Fig. 3. — Expansion ratio of half-mass radius as a function of SFE. (A color version of this 
figure is available in the online journal.) 
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Fig. 4. — Bound mass fraction as a function of SFE. (A color version of this figure is available 
in the online journal.) 
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3.2. Cluster and cloud have different initial density profiles 

We suppose that it is too idealized to consider both cluster and cloud have the same 
initial density profile (as in the previous section), and it is conceivable that they are more 
likely to have different initial profiles. To study the cases with different initial profiles, we 
fixed the Plummer length scale of the cluster at 0.6 pc, but varied the Plummer length scale 
of the cloud from four times smaller to four times larger than that of the cluster. Thus, there 
are two parameters to describe the initial profile of the cloud: {Mf,,ao}, its mass and the 
initial Plummer length scale. (One can use the SFE t] and oq instead.) It is obvious that rj 
alone is not sufficient to describe the results with different Oq and Mf, (or equivalently rj, see 
Equation ([T])). When we plot the expansion ratio e of 20 Oq cases against r] (see Figure E]), 
no discernible relation(s) can be found. (We should point out that if we plot the result of 
one Co, the graph will be similar to Figure O) 

However, e does show a "simple relation" with {Mh,ao}. Figure [H] shows the contours 
of e in the parameter space {Mb, qq} for the dispersing timescale te = 3.3 Myr. Under- 
standably, e is large/small (cluster is destroyed/remains intact) when oq is small/large and 
Mb is large/small (or t] is small/large). In fact, there are three final configurations for the 
cluster: (i) destroyed, (ii) loose structure, and (iii) compact core. The boundaries between 
these configurations closely mat ch the contours e = 10 and 2 for the cases with different 
dispersing rates (for details, see IChen fc Koll2008h . In other words, on the one hand for a 



fixed effective size of the natal cloud, it is clear that the survivability of the cluster is higher 
when the SFE is larger. On the other hand, for a fixed SFE, the survivability increases when 
the effective size of the natal cloud increases. For example, when the cloud contains 90% 
of the mass of the system (i.e., total SFE 10%), the half- mass radius of the cluster expands 
at least twice in 30 Myr. The smaller the size of natal cloud, the more the star cluster will 
expand. Apparently, the reason is the dynamics of the cluster is dominated more by the 
cloud when the mass of cloud initially enclosed within the cluster is more, such as the cases 
of small SFE and small effective natal cloud size. And in contrast for the cases of large SFE 
and large effective natal cloud size, the dynamics is dominated by the cluster. 

Although SFE t] fails to be a good survivability indicator of the two-parameter system, 
we are able to identify other survivability indicators: the initial cluster-cloud mass ratio and 
the initial kinetic energy. We describe these two indicators in the following separately. 

Cluster- cloud mass ratio 

First denote the mass fraction / Lagrangian radius of the cluster as rj, i.e., within rj the 
mass of the cluster is fM^, and denote this mass by M^^vf = fM^. Also denote Mb^rf as the 
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Fig. 5. — A plot of the expansion ratio of half-mass radius against SFE for 20 cases of oq. 
(A color version of this figure is available in the online journal.) 




Fig. 6. — Contours of the expansion ratio of half-mass radius (erft,„ = Crso? solid line) in the 
parameter space {Mf,,ao}. The dashed lines are contours of the cluster-cloud mass ratio at 
half-mass radius, Pr^^ = Prso- ^^i^ example, the cloud dispersing timescale is te = 3.3 
Myr. The corresponding SFE is labeled at the top. (A color version of this figure is available 
in the online journal.) 
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mass of the cloud within Vf initially. Define the cluster-cloud mass ratio as (see lChen Sz Ko 



2008|) 

This closely resembles the definition of SFE t] (see Equation ([T])), but Prj. contains information 
about the length scale of the cloud Qq- Note that Poo = V (when / = 1, — > oo). Since we 
are using the Plummer model for both cluster and cloud, there is a simple relation between 
Prf and 7] 



where / = M^^rf/Mc and F = Mb^rf/Mf, = ry{rj + alY^"^. We also extend the half-mass 
radius expansion ratio of Equation (jlj) to 

_ r/ (final) 

r^initial)- ^' 



We find that for each dispersing timescale t^, f can be tuned such that data from the 
two-dimentional parameter space {M^, oq} collapses to a relation in the f3rf — plane, see 
Figure [7l The relation closely resembles the relation rj — e for cases of fixed oq (compared 
to Figures [7] and [5l) and the turnover points for te = 3.3,1.1,0.65,0.33 Myr are at jSrj ~ 
0.18,0.28,0.40,0.53, respectively. We note that r/ does not increase monotonically with 
increasing dispersion timescale. 

We stress that although it may seem a bit contrived that we have to tune rj to obtain a 
tight relation, the result is reasonably good even if we pick one particular rj for all timescales. 
Figure E] shows the relations at half-mass radius for all dispersing timescales ( gainst 
/3r^^). In Figure [6l we see that the contours of e,.^^ and Pr^^ agree with each other in the 
parameter space {Mh,ao}. 



Initial kinetic energy or virial energy 

To clarify the relation between the final state and the initial condition, we also check the 
energies. We find that the initial kinetic energy is nicely correlated with the expansion ratio 
of half-mass radius e, see Figure |9l Since the initial system is in virial equilibrium, the initial 
virial energy is just twice of the kinetic energy, hence it should also be a good indicator. 

The initial kinetic energy and the initial virial energy show how the cloud mass affect 
stars inside the cluster. is the cluster-cloud mass ratio within a certain region. All of 
them show nice relations with the expansion ratio e or e,, . But if one would like to apply to 
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Fig. 7. — Expansion ratio of the chosen Lagrangian radii as a function of the cluster-cloud 
mass, (3rf (A color version of this figure is available in the online journal.) 
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Fig. 8. — Expansion ratio of the half-mass radius as a function of the cluster-cloud mass, 
(3rf Here / is 0.5. (A color version of this figure is available in the online journal.) 
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observations, it might not be easy to obtain the energy estimation for the whole cluster. (3rf 
should be the more convenient choice. 

Figure [TU] shows three initial mass density profiles of clusters and clouds at a particular 
cloud mass M^/Mc or total SFE t]. In the figure, M^/Mc = 4 and r] = 0.2. The upper 
panels show the density profiles against radius and the lower panels show the corresponding 
(5r against radius. Note that j3r —>■ t] as r ^ oo. In the middle column, = clc, the two 
density distributions are the same, and the corresponding is always 0.2 at any radius. In 
the left column, ao < ac (i.e., comparatively the cloud mass is distributed more concentrated 
at the origin than the cluster), i3r increases monotonically with r (from < rj towards r] as 
r ^ oo and l3r always smaller than rj). In the right column, oq > Cc (i.e., comparatively the 
cloud mass is distributed more evenly than the cluster), jSr decreases monotonically with r 
(from j3r > rj towards 77 as r — oo and (3r always larger than 77). Therefore, for every total 
SFE, we have three distinct types: oq <, =, > ac- The cluster has a higher/lower chance 
to survive when the cloud size is larger /smaller than the cluster size at a given total SFE. 
This also suggested that at a low total SFE the cluster has a higher chance to survive if the 
length scale of cloud is larger. 



To this end, we would like to point out a similar idea introduced by iGoodwid (120081 ). 
The effective SFE, eSFE, is defined as the virial ratio of the stars after instantaneous gas 
expulsion. The total SFE equals the eSFE when the cluster and cloud share the same 
density profile. When the cluster is born dynamically "cold/hot", the velocity of a cluster 
is smaller /larger than what is required for virial equilibrium, the eSFE is larger /smaller and 
the cluster has a better/worse chance to survive. This is comparable to our results of clouds 
with larger /smaller length scales. 

Moreover, eSFE is closely related to the initial kinetic energy of the cluster that we 
mentioned earlier. For a fixed total SFE, the initial kinetic energy is smaller if the size of 
the natal cloud is larger. In this case the cluster is 'colder' and tends to survive after the 
cloud dispersed. 



3.3. Mass segregation 

Mass segregation is another important issue in stellar clusters. Once again we use 
NB0DY2 to study mass segregation. All simulations start with the Salpeter mass function, 
i.e., the mass function slope is —1.35. We find that /J^^^ is also a good indicator for the level 
of mass segregation. Figure [11] shows the mass function slope derived from within the final 
half density radius as a function of [3r^^ for cases with e^.^^ < 20. When [3r^^ is small (i.e.. 
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Fig. 9. — Initial kinetic energy and the expansion ratio of half-mass radius. The units of the 
initial kinetic energy here are in iV — body units. (A color version of this figure is available 
in the online journal.) 



cloud mass is the dominant factor), the mass segregation of the cluster is less severe. We 
note that using different f^rf does not change the relation much. 



The preferential loss of low- mass stars is due to two-body relaxation (see ISpitzerlll958 



Spitzer fc Harmlll958l ) and would flatten the mass function. We should point out that we fit 
the mass function by a simple single power law, which is not the be st fit for all case s. As a 
matter of fact, the broken power law observed in the Arches cluster (IKim et al.ll2006l ) is also 
seen in our simulations. We no te that this could also be ob tained by dynamical evolution 
without a natal cloud (see also iPortegies Zwart et al.l 120071 ). However, we cannot find an 
obvious relation between the broken power law and SFE t] or the cluster-cloud mass ratio 
j3. This may be due to poor statistics in our runs (2500 stars initially). Moreover, the broad 
distribution of the final mass function slope in Figure [11] is perhaps a result of the simple 
single power law fitting. 



4. Summary and Remarks 

We studied the consequence for clusters after parent clouds disperse in different dispers- 
ing timescale, cloud mass, initial cloud length scale, and IMF. We mimic cloud dispersal by 
increasing cloud size or length scale, but keeping the total cloud mass unchanged. 



-14- 




Fig. 10. — Mass density profiles and the corresponding cluster-cloud mass ratio, for cases 
with different length scales of the cloud. The total SFE 77 = 0.2 (i.e., M^/Mc = 4). From left 
to right, the ratios of length scale of the cloud to cluster are 0.5, 1, and 2, respectively. The 
mass of cloud is more concentrated (cloud-dominated) in the left column and more spread 
out (cluster-dominated) in the right column. The mass densities are in arbitrary units. 
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Fig. 11. — Mass function slope as a function of half-mass (3 value. The mass function 
slope is derived from stars within the final half-density (2D) radius. The IMF slope is —1.35 
(Salpeter) for all cases. (A color version of this figure is available in the online journal.) 

SFE, r], is a good indicator for the survivability of embedded cluster after the dispersal 
of the parent cloud only if the cluster and cloud have the same density profile. It fails badly 
if the cluster and cloud have different density profiles. The latter cases can be described 
by two parameters {M^, oq}, the mass of the cloud and the length scale of the Plummer 
model. The expansion ratio of the cluster is obviously depending on the two parameters. 
Although we cannot use 77 as a survival indicator, we find two quantities that can serve the 
same purpose. 

The two quantities are (i) the initial cluster-cloud mass ratio (within a certain La- 
grangian radius of the cluster), (see Equation ((51)), and (ii) the initial kinetic energy (or 
virial energy). Each of them can organize the data of itself and the expansion ratio (of a 
certain Lagrangian radius of the cluster) at a dispersing timescale into a 'one-dimensional' 
relation (effectively collapsing data from a surface into a line). Moreover, the relation is 
very nice in the sense that it has two branches, one corresponds to cluster destruction and 
the other to survival. Hence the two quantities are good indicators for the survivability of 
embedded clusters. 

We should point out that for different cloud dispersing timescales we need to choose 
somewhat different Lagrangian radius to obtain the tightest correlation between the cluster- 
cloud mass ratio and the expansion ratio. Having said that, in fact, P,.^^ fits the data rea- 
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sonably well for all different cloud dispersion rates . The tight correlation strongly suggested 
that the evolution of the embedded cluster depends on the total mass enclosed within the 
corresponding Lagrangian radius of the cluster. We thus propose that the relevant (initial) 
effective timescales should be defined in terms of the total mass enclose d within the best-fit 
Lagr angian radius (or roughly the half- mass radius where / = 0.5) (e.g., Binney &: Tremainel 
1987h . 



''"relax 



fN 



where is the number of stars in the cluster, and Trciax and Tcross are, respectively, the 
effective relaxation time and crossing time. 

Furthermore, for clusters with Salpeter mass function, the slope of the final mass func- 
tion shows a roughly linear relation with the cluster-cloud mass ratio at half-mass radius, 
Pr^^, see Figure [TTJ In addition, all the dispersing rates we considered share the same 
relation. 

As mentioned in Section [21 similar works have been done with an e qual-mass model 



[e.g.. lGeyer fc Burkertll200ll : iBaumgardt fc Kroupall2007l ). iGoodwinI (119971 ) mentioned that 
the survivability of the clusters is only slightly higher in the cases of equal mass than cases 
with mass function. For comparison with the result of mass function presented here, we also 
did the same simulation survey for the cases of equal mass. We confirm that there is no 
significant difference on the total mass of escaping stars. However, the mass distributions 
would be different for radii larger than ryo (70% Lagrangian radius) in these two models. 

Observationally, the initial kinetic energy of the whole star cluster is not easy to estimate. 
On the other hand, the cluster-cloud mass ratio is more promising and should warrant further 
studies. 



We thank the referee for valuable suggestions to an earlier version of the paper. This 
work was supported in part by the National Science Council, Taiwan under the grants NSC- 
96-2112-M-008-014-MY3. 
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